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Abstract. We introduce a spectral transform for the finite relativistice Toda lattice 
■ (RTL) in generalized form. In the nonrelativistic case, Moser constructed a spectral 

pH , transform from the spectral theory of symmetric Jacobi matrices. Here we use a non- 

symmetric generalized eigenvalue problem for a pair of bidiagonal matrices (L, M) to 
define the spectral transform for the RTL. The inverse spectral transform is described in 
terms of a terminating T-fraction. The generalized eigenvalues are constants of motion 
and the auxiliary spectral data have explicit time evolution. Using the connection with 
the theory of Laurent orthogonal polynomials, we study the long-time behaviour of the 
\jr^ | RTL. As in the case of the Toda lattice the matrix entries have asymptotic limits. We 

lO ■ show that L tends to an upper Hessenberg matrix with the generalized eigenvalues 

sorted on the diagonal, while M tends to the identity matrix. 

o 

CI 
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^ ■ 1. Introduction 

The relativistic Toda lattice (RTL) was introduced by Ruijsenaars |l7j and studied in 
0, ||, f|, [14], |18[ , for a review see [|12|] . The finite RTL is defined by the system of equations 



c3 



> 



X ' •• _ 2 • ( ■ exp(g n _i - q n ) . exp(g ra - q n+l ) 



n V " _1 l + e 2 exp(g n _i - q n ) n+1 1 + e 2 exp(q n - q n+1 ) 
with iV 6 N, 1 < n < N and the convention q = — oo and q^+i = +oo. Here the 
q n are functions in the time parameter t. The system ( |1 . 1|) is the Newtonian form of a 
Hamiltonian system with Hamiltonian 

N 

H(q 1 , ...,q N ,p u ...,p N ) = J^e Pn /i(g n _i - q n )h(q n - q n+1 ), 

n=l 



where h(x) = vl + e 2 e x . Setting q n = Q n + 1/e and letting e — 0, one can easily check 
that ( |1 . 1| ) reduces to the equations of motion of the finitie nonrelativistic Toda lattice, 

Q n = exp(Q n _! - Q n ) - exp(Q n - Q n+1 ), 1 < n < N, 

so that ( |1 . 1| ) is a one-parameter perturbation. 
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In analogy with the Flaschka variables for the nonrelativistic Toda lattice one can 
use a (non-invertible) change of variables to prove integrability. Bruschi and Ragnisco 
|], |2|, see also |12|, |TBf , obtain the two forms 

• _ fe n _ b n-l 



O-n+1 O-n-1 



b n i 



On O-n+1 



1 < n < N, 
1< n < N - 1, 



[1-2) 



and 



O-nibn-1 — On 



1 <n < N, 
1< n < N- 1, 



(1.3) 

= b n (a n — a n+ \ + a n _i — a n+ i), 

both with a n > for 1 < n < N, b n > for 1 < n < N — 1, and b = 0, b N = 0. 
The systems (|1.2j) and (|1.3| ) can also be written in matrix form. Define two bidiagonal 
matrices L and M by 
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a 3 ' 
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Then Suris |T8|1 noted that fll 
L 



can be written in the Lax form 



LA - BL 
M = MA- BM, 



1.5) 



where 
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1.6) 



Then A = — (L X M)_ and B = —(ML x )_, where we use X_ to denote the strictly 
lower triangular part of X. The other system ( |1.3| ) can also be written in the form ( |1 . 5| ) 
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but now A = 
we show that 
namely 



see Remark 13.41 below. In Section 3.1 



(M~ 1 L)_ and B = —(LM~ r 
172) and (|1.3|) are special cases of a generalized form of the finite RTL, 



L = 
M 



F (LM _1 ) 
F (LM' 1 



l-l(f {m- x l) 

M - m(f(M- 1 L) 

. For F(x) 



(1.7) 

1/x, we obtain ( |1.2| ) 



where F : R + — > R is an arbitrary function on IE 
and for F(x) = x, we obtain (|1.3j ). 

The main goal of this paper is to solve the finite RTL in its general form (|1.7|) with 
the aid of a direct and inverse spectral transform, in a similar way as was done by Moser 
13| for the nonrelativistic case. While the spectral transform for the nonrelativistic 



case is based on the spectral theory of tridiagonal (Jacobi) matrices and is connected 
with orthogonal polynomials, we will show that in the relativistic case, the spectral 
transform is based on the spectral theory of pairs of bidiagonal matrices ( |1.4| ) and 
Laurent orthogonal polynomials jR], [16], [2D], ||. See Figure [I] for the general scheme of 
the direct and inverse spectral problem. 



matrix data 
at t = 



matrix data 
at t 



direct 
spectral 
problem 



inverse 

spectral 

problem 



spectral data 
at t = 



evolution of the 
spectral data 



spectral data 
at t 



Figure 1. The scheme of the direct and inverse spectral problem. 

To the best of our knowledge, the only transform methods that have been explicitly 
introduced in the study of the RTL are the ones by Ragnisco and Bruschi JITJJ for the 
periodic RTL, by Ruijsenaars [17], and the scattering transform of Cosentino || for the 
infinite RTL. 

The Laurent orthogonal polynomials P n are defined by 

P„ = l, P n (z)=det(zM n -L n ), l<n<N, (1.8) 

where L n and M n are the n x n upper left corner blocks of L and M, respectively. They 
satisfy the recurrence relation 

P n (z) = (z - a n )P n ^i(z) - b n -\zP n -i{z), 1 < n < N, 
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with P = 1, P_i = 0. Under the assumption that a n > and b n > 0, it is known 



that the zeros of these polynomials are simple, real and positive, see |9|, [K| [11], [L6| and 
Lemma [2.1| below. The first set of spectral data we choose are the ordered zeros of the 
polynomial Pn(z) = det(zM — L) which we denote by < Ai < ■ • • < Ajy. Then the Xj 
are also the eigenvalues for the generalized eigenvalue problem of the pair (L, M) 

Lx = XMx. (1.9) 

To obtain our second set of spectral data we associate with each A, the value 



w i 



(ujMvi)' 



where uj and Vj are the left and right eigenvectors of the generalized eigenvalue problem 
( |1.9|) , respectively, normalized to have their first component equal to 1. We prove in 



Section ^J] that Wj > and J2f=i w j = 1- The direct spectral transform 

ai,...,ajv>0 1 j < Ai < ■ • ■ < Aat 

&i, . . . , 6jv_! > J ) w h . . . , w N > 0, YlJ=i w j = 1 



is invertible. Using T-fractions || [TT| , we explain in Section ^T|] how to obtain the a n 
and b n when the spectral data Aj and Wj are given. We also note that the polynomials 
P n satisfy 

N 

V P n (A i )P w (A i )^ = 0, if n > m, 
j 



which is known as Laurent orthogonality P, |^, [10|, [16] 



In Section p.2| we describe the time evolution of the spectral data for the finite RTL 
in the generalized form (|1.7f) with positive initial data. We show that the Aj are time 
independent and that 

w iW = -tt^ , 1<J<N. 

k=l 

Finally, in Section § we investigate the behaviour of the RTL for t — > ±oo. We 
assume that F is strictly increasing. Then 

lim b n (t) = (1.10) 

t— >±oo 

and 

lim a n (t) = Ajv-n+ij b m a, n (t) = A n , 1 < n < N. (1-H) 

Similar limits hold for the case that F is strictly decreasing. This sorting property of 
the RTL indicates that the choice of the spectral data A x < • ■ ■ < Ajv is very natural. In 
Theorems L3 and L4 we establish the precise rates of convergence of (|1.10|) and ( 1 . 1 1|) . 



This is similar to Deift et al. || who obtained precise rates of convergence for the SVD 
flow (Singular Value Decomposition). 
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2. The direct and inverse spectral problem 

In this section we introduce a transformation from matrix data a n > 0, b n > to spectral 
data Xj, Wj. We prove that this spectral transform is a bijection between (M + ) 27V_1 and 
its image {0 < Ai < • • • < Ajv, Wj > 0, J2f=i w j = !}■ 



2.1. The direct spectral problem 

We start from positive matrix data a n , 1 < n < N, and b n , 1 < n < N — 1, and 
construct the bidiagonal matrices L and M as in ( |1.4| ). 
Define the finite set of monic polynomials P n as 

P = 1, P n (z) = det(zM n - L n ), l<n<N 

where L n and M n are obtained from the matrices L and M, respectively, by deleting 
the last N — n rows and columns. These polynomials satisfy the recurrence relation 

P n (z) = {z- a n )P n ^(z) - b n ^zP n _ 2 (z), l<n<N, (2.1) 

as can be easily seen by expanding the determinant det(zM n — L n ) first by the last 
column and then by the last row. The initial conditions for (|2.1| ) are P Q = 1 and 
P_i = 0. The zeros of these polynomials behave in a way which is similar to the 
behaviour of zeros of polynomials which are orthogonal on the positive real line. The 
following lemma can be found in P, |10|, O, [TH], but we include a proof for completeness. 



Lemma 2.1 The zeros of the polynomials P n (z) = det(,zM„ — L n ) are real, simple, 
positive (0 < z\ )Tl < ■ ■ ■ < z n , n ) and have the interlacing property, which means that 

< 

Proof. The lemma is true for the case n = 1 because Pi{z) = z — a± and a\ > 0. 
Suppose that the lemma is true for some n < N. By the interlacing property we then 
know that P n -i changes sign between any two consecutive zeros of P n . Evaluating the 
recurrence 

P n +i{z) = {z- a n+ i)P n {z) - b n zP n _i(z) 

at a zero z j>n of P n , we get P n+X (zj >n ) = -b n z^ n P n -\{z i{n ). Since b n > and Zj >n > 0, we 
see that P n+ \ and P n _i have opposite signs at the zeros of P n . It follows that also P n+ \ 
changes sign between any two consecutive zeros of P n , and therefore it has at least one 
zero in each of the intervals (zj in , Zj + i jn ), 1 < j < n — 1. Since P n -i(z n ,n) is positive, it 
also follows that P n+ i(z n ^ n ) is negative. Since P„+i is a monic polynomial, it must have 
a zero in (z n>n , oo). 

To complete the proof we will show that there is also a zero in (0, Zi iTl ). Evaluating 
the recurrence at 0, we get P n+ i(0) = — a n+ iP n (0). Since P n is a monic polynomial of 
degree n having all its zeros on the positive real line, we find that (— 1)"P„,(0) > 0. Thus 
(— l) n P n+1 (0) < 0. Since P„_i is a monic polynomial of degree n — 1 with all its zeros 
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to the right of z ln , we obtain (— l) n P n -i(zi tn ) < 0. Since P n +i and P n -\ differ in sign 
at zx >n , we get (— l) n P n+ i(,2i in ) > 0. So P n+ \ has opposite signs in and Zi >n , which 
shows that there is indeed a zero in (0,2;i jn ). The lemma now follows by mathematical 
induction. □ 



Remark 2.2 An analogous reasoning as in Lemma |2.1| can be applied to the 
polynomials A n (z) = det(zM n — L n ), A = 1, where L n and M n are obtained from 
the matrices L and M by deleting the first N — n rows and columns, respectively. These 
polynomials satisfy the recurrence 

A n (z) = {z- aAr_ n+ i)A n _i(» - b N ^ n+1 zA n ^ 2 (z), (2.2) 

for 1 < n < N, with initial conditions Ao = 1 and A_i = 0. So also the zeros of these 
polynomials are real, simple and positive and have the interlacing property. 

Consider now the generalized eigenvalue problem of the pair (L, M) 

Lx = XMx. (2.3) 
The eigenvalues are the zeros of the polynomial P/v(A) = det(AM — L). From Lemma 



27T| we know that these eigenvalues are real, simple and positive. We denote them by 

< Ai < • • ■ < X N < oo. 

We use Vj to denote the right eigenvector corresponding to the eigenvalue Xj , normalized 
so that the first component is equal to 1. From the bidiagonal structure of the matrices 
M and L it easily follows that the first component of any eigenvector of ( |2.3| ) is non-zero, 
so that this normalization is always possible. We use Uj to denote the corresponding 
left eigenvector, also normalized to have the first component equal to 1. We view Uj as 
a column vector so that ujL = XjtlJ M. We collect the eigenvectors in the matrices 

V = (vf - - vn) and U = {vi\ ■ ■ ■ un) ■ 

A left and right eigenvector of this generalized eigenvalue problem, corresponding to 
different eigenvalues, are M-orthogonal, that is, uj Mvu = whenever j ^ k. This 
implies that U T MV is a diagonal matrix. Since the three matrices U, M, and V are 
invertible, the diagonal elements of U T MV are non-zero, and we can make the following 
definition. 



Definition 2.3 Let U and V be as above. Then we define 
(U T MV)jj 

Now we have defined the spectral data Aj and Wj for 1 < j < N. To complete the 

iat in. : "> fl anH V . 

■T- 



^ = 777t^7Y-=(^) > 1<J<N. (2.4) 

Xj and Wj for 1 < j < N. 1 
description of the direct spectral transform we show that Wj > and Yl!j=i w j = 1- To 
prove this we introduce the Weyl (or resolvent) function 

f(z) = e?(zM - L)~ l e\ = [{zM -L)' 1 )^, zeC\ {A 1; . . . , A^}, 
with e*i the first unit vector in ~R N . 
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Lemma 2.4 The Weyl function is given by 

3=1 J 

Proof. Note that by definition of Wj, we have uj Mv^ = —S-jk, and so 
ujLv k = XjujMvk = —S jjk , 1 < j, k < N. 

Wj 

In matrix form this is U T MV = diagf — —) and U T LV = diagf . . . , 
Combining these two equations, we find that 

U T (zM-L)V = diag( Z -^,..., Z -^). 

\ Wi w N / 

From this we readily obtain that 

(zM - Ly 1 = Vdiagf^-, . . . , -^-)u T . (2.5) 

\z — Ai z — Aat/ 

We now recall that the eigenvectors Uj and Vj are normalized with their first component 
equal to 1. Thus the first rows of the matrices U and V contain only ones. Then ( |2.5| ) 
gives that 

N 

/(,) = ((,*-*)-') =£-SL. 

j = l J 

□ 

The lemma shows that the Xj are simple poles of the Weyl function and that the 
Wj are the corresponding residues. A second representation of / is in terms of the 
polynomials A n introduced in Remark [2.2| . 

Lemma 2.5 The Weyl function f can be written as f(z) = ^ 1 ^ . 

Proof. This follows immediately from Cramer's rule. □ 

As a corollary we get 
Corollary 2.6 We have Wj > for every j and Y^,j=i w j — 1- 
Proof. If we combine Lemma [2.4| and Lemma |2.5|, we get 



^N-l(z) = y> Wj 



From Remark |2.2| we know that the zeros of A^-i and A^v are interlacing. Since Ajv-i 
and An are monic polynomials, it then follows that Wj > for every j. The fact that 
J2j=i w j — 1 follows by multiplying ( ^6|) by z and then letting z — > oo. □ 
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We now have the direct spectral transform 

(R+) 2JV 1 ->• A: (oi,...,ajv,6i,...,6jv-i) >-> (Ai, . . . , X N , w x , ■ ■ ■ , w N ) 

where 

< Ai < • • ■ < X N < oo, 



A = < (Ai, ...,X N ,wi,... ,w N ) 



w i > °> Eli w i = 1 



(2.7) 



We conclude this paragraph by establishing an explicit representation of the Weyl 
function in terms of the given matrix data. This representation is also the key to the 
inverse spectral transform. It turns out that the function zf(z) can be written as a 



special continued fraction, known as a T-fraction, see pi 
Lemma 2.7 The Weyl function f can be written as 

/(*) = ^# = — h • (2-8) 

b 3 z 



z - a 2 - 



z - a 3 - 

b N z 



z — a A? 



Proof. Let r k (z) = ior 1 < k < N. Since A (^) = 1 and Ai(z) = z — we 

have ri(z) = z _ a . From the recurrence (|2.2|) satisfied by the polynomials A^, we get 
that 

r k( z ) = ^7 7 7- 

Z — OjV-fe+1 — Z0N-k+2 r k-l{Z) 

Then it is easy to show by induction that 

v ' z — dN-k+i r ; 

z — GtjV-fc+2 



z — a N _ k+3 - 



b N z 



z — a, 7v 

Taking k = N in ( [2.9|) we obtain the lemma. □ 



From Lemma p.4| and Lemma \L7\ it is easy to see that the spectral transform is 
injective. In the next section we show that it is also surjective. 



2.2. The inverse spectral problem 

We now start from spectral data (Ai, . . . , Atv, Wi, . . . , wtv) in A, see ([2.71) , and show that 
there exist positive matrix data (ai, . . . , a/v, b\, . . . , 6jv-i) which correspond to the given 
spectral data. We also give a method to construct these matrix data. In view of Lemma 
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[2.7| it is our task to develop Ylj=i J^\~ a continued fraction of the form ( |2.8| ) with 
positive a n and b n , that is, we want 

%^r— l — • (2 - 10) 

7=1 J Z — Cli — r 

3 b 2 z 

z - a 2 



z - a 3 - 

bN-lZ 



Z — <2jv 



Theorem 2.8 For given spectral data (Ai, . . . , Ajv, w\, . . . , wn) in A there exist positive 
a n and h n such that ( |2.10| ) holds. 



Proof. Let us say that / G Wn if / is of the form 



N 

j=i j 



with < Ai < • • • < \n, Wj > 0, and J2f =1 Wj = 1. We will prove that for / G Wjv, 
j >( 

m 



'3 

there exist a,b > and g G Wn-i so that 

1 



2 — a — bzg(z) 

The theorem then follows by repeated application of this fact. 

To prove the claim, we observe that / is a rational function which we can write as 

with Pn-i arid monic polynomials of degrees N — 1 and N, respectively (here we use 
the fact that J2f=i w j = !)■ Then there exist a, b G R and a monic polynomial pn-2 of 
degree at most N — 2 such that 

gjvO) = (z - a)p N -i(z) - bzp N ^ 2 (z). (2.12) 

Taking z = in ( p.!2|) we find 

9jv(0) / ^ wj 



N ^ - 1 



a 



Pjv-i(O) \^ 



which is positive, since all Xj and Wj are positive. To find b, we write ( |2. 12| ) as 



—q N (z) + [z — o)pn-i(z) = bzp N _ 2 (z) and compute the coefficient of z N 1 on the 
left-hand side. Since qis{z) = Yl^ = i(z - \j) and pn-i{z) = J2k=i w k Hj^k( z ~ -\/), tnis 
coefficient is 

N N N N / N \ 

A ^ ~ u,fc a j ~ a = X A 1 ~ Wk )- a 

3=1 k=l j=l,j^k j=l \ fc=Wj / 

N / JV x _1 

10 1 
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In the last step we used ( |2.13| ) and Y2f=i w j = 1- Next, from Jensen's inequality 
zdfM(z) I < / (p(z)dfi(z) 



applied to the strictly convex function <j>{z) = l/z and the probability measure 



N / N \ -1 

i=i \i=i V 

So the coefficient of z^ -1 in bzpN-2{z) is positive. Since ptv_2 is a monic polynomial of 
degree < N — 2, it then follows that pn-2 has exact degree N — 2 and that b is equal to 
Q2.14Q . Hence b is positive. If we now define g = PN-2/PN-1 then we have that 

f(z) = ; — -, — r, with a, b > 0. 

z — a — bzg(z) 

It remains to show that g G Wjv-i- Observe that the poles of (7 are the zeros of 
/. From ( p. 11] ) with Wj > it follows that the zeros of / interlace with its poles 
< Ai < A2 < • • • < \n- So the poles of g are simple and positive. If we denote them 
by < A* < A£ < ■ • • < A^_ 1( then 

/ \ N-l * 

p N _ 2 (z) 



i=i ? 

AT-l 



for certain residues u>* that satisfy Ylj=i w* = 1, because Pn-2 and pat_i are monic. 
We also have 

,,, , , , z — X) f z — a 1 \ 
Wj = tan (, - X M z) = tan — ( — - j 



z-A* -l 



lim 



z^bzf(z) bxmxy 



N 



Since b > 0, A* > 0, and /'(A*) = — 77 t~t^ < 0, we see that w* > for every j. 

This completes the proof of the theorem. □ 



2.3. Laurent orthogonality 

Suppose that \l is a positive measure with support in M + , for which all the strong 
moments J z k dfi(z), fcGZ, exist. For n not exceeding the number of mass points in the 
support of /i, the monic Laurent orthogonal polynomial, p n , of degree n is defined as the 
monic orthogonal polynomial of degree n for the varying measure dfi n (z) = z~ n dfi(z). 
So p n satisfies the orthogonality relations 

j Pn ( z )z k ^l = 0, 0<*<n-l. 
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Cochran and Cooper J| and Jones et al. || ||, |10| showed that these polynomials satisfy 
a recurrence relation of the form 

p n (z) = {z- a n )p n _i(» - I3 n _izp n _ 2 {z), n > 1, 

with initial conditions p = 1 and p_\ = 0, and positive recurrence coefficients a n and 
f3 n -\. From the theory of orthogonal polynomials it follows that the zeros of p n are 
simple, real and positive, since fi is a measure on the positive real line. Some examples 
of (continuous) Laurent orthogonal polynomials are given in 

We now show that the polynomials P n (z) = det(;zM n — L n ) introduced in (|L~8| ) 
are the monic Laurent orthogonal polynomials with respect to the discrete probability 
measure 

N 
3=1 



see also |16| for a similar situation. These polynomials are also generated by the 
recurrence relation (|2.1|) . This illustrates the connection between the spectral transform, 
proposed in Section |2.1], and the theory of Laurent orthogonal polynomials. 



Theorem 2.9 Let a n > 0, 1 < n < N, and b n > 0, 1 < n < N — 1. Let the polynomials 
P„ be defined by the recurrence relation fl2.1| ) with P = 1 and P_i = 0. Let Xj and Wj 
be the spectral data associated with a n , b n as defined in Section |2"X Then 

N n 

J]P m (A j )P n (A J )^ = 5 mtn l[b k , 0<m<n<N -1. 

j=l i k=l 

Proof. Evaluating the recurrence relation ( |2.1| ) at Xj (which is a zero of P/v) one can 
easily see that the right eigenvector Vj for the generalized eigenvalue problem ( |2.3| ) is 
equal to 



V j =(l,F L (\ j ),...,P S - 1 (\ j )) . (2.15) 
Also note that the left eigenvector Uj can be written as 

/ Pi(Aj) P 2 (A,) Pv-i(Aj) V 

3 ~ {' Xjbi ' Xp 2 bi ' • • • ' Af-^-i ...h) ' 1 ' 

Since U T MV = diag^^-, . . . , ^-j by the definition of the Wj, we have 

I/diag(wi, . . . , w N )U T = M _1 . 

Now note that M is lower triangular with ones on the diagonal, and therefore the same 
holds for M _1 . SoifO<m<n<A^— lwe have 

( V^diag(w;i, . . . , w N )U T ) = 8 mn . 

V / m+l,n+l 

Using ( |2.15| ) and (|2.16|) we finally get 

N 

^ ^ P?n ( A j ) Pi ( A j ) — Q~m,ni 

i=i A" n ^ 

fc=l 
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which proves the theorem. □ 

The theorem suggests an alternative method for performing the inverse spectral 
transform. Knowing the spectral data, one can construct the monic Laurent orthogonal 
polynomials corresponding to the discrete measure /ijv- These polynomials satisfy a 
three term recurrence relation with coefficients equal to a n and b n . 

3. The generalized finite RTL and spectral evolution 

In the introduction we have already mentioned that the finite RTL 
a n = - 1< n < N, 

" a-n+l a n -x ' — — ' 

b n = b n ( i - , 1< n < N - 1, 
with 6q = 0, 6at = 0, can be written in Lax form representation Hl8| , namely 



L = LA — BL 

(3.1) 

M = MA- BM, 

with A = — (L -1 M)_ and B = — (A/L _1 )_. Here the matrices L and M contain the 
matrix data a n and b n and are defined in ( |1.4| ). In this section we first show that this 
system can be generalized as in (1.7) . Then we will solve the finite RTL in its generalized 



form, using the spectral transform introduced in Section |^. We will find the explicit 
evolution of the spectral data. 

3.1. The generalized finite RTL 

The generalized finite RTL is defined by the matrix differential equations 
L—(f (LM" 1 )) L-l(f {M~ l L) 



M=(F (LM" 1 )) M-M^F (M~ l L) 



(3.2) 



where F : IR + — > M is a real- valued function on IR + . The matrices L and M contain the 
matrix data a n and b n as in (|1.4j) . We assume positive initial data a„(0) > 0, b n (0) > 0. 
The system is of the form (|3.1| ) if we set 

A = -(f (M _1 L)) and B = -(f (LAr 1 )) . (3.3) 

The eigenvalues of LM~ l and M~ X L are equal to the generalized eigenvalues for 
the pair (L,M). We know that there are N distinct positive generalized eigenvalues. 
Hence the matrices LM~ X and M~ l L are diagonalizable, and in fact we have 

LM' 1 = U- T DU T , M~ l L = V^DV 

where D is a diagonal matrix containing the generalized eigenvalues Ai < A2 < ■ • • < Xn- 
Thus 



F (LM' 1 ) = f/- T diag(F(A 1 ), . . . , F(X N ]^U 1 
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and 

F (M^L) = V-Miag^Ax), . . . , F(X N fjV, 

see 0, Definition 6.2.4]. 

Of course one has to check that the structure of the matrices L and M is preserved 
by the differential equations ( |3.2| ). We will use the following lemma. 

Lemma 3.1 For every function F we have that 

LF (M _1 L) = F (LM- 1 ) L and MF (M~ l L) = F (LM' 1 ) M. (3.4) 

Proof. For a monomial F(x) = x n , we have 

LF (M _1 L) = L (M~ l L) n = (LM- 1 ) 11 L = F (LM' 1 ) L 

and 

MF (M~ X L) = M (M _1 L) n = (LM~ l ) n M = F (LM' 1 ) M. 

Then, by linearity, ( |3.4|) holds for every polynomial F. On the finite set of eigenvalues 
of LM~ l (which are the same as the eigenvalues of M~ l L) each function F is equal to 
a polynomial, so the equalities hold for every F. □ 

We now prove that the proposed generalisation of the finite RTL is justified. 

Theorem 3.2 The system of matrix differential equations ( |3.2| ) is well defined, which 
means that LA — BL is a diagonal matrix and MA — BM has only non-zero elements 
on the first subdiagonal, where 

A = -(F (M -1 !,)) and B = -(f (LM" 1 )) . 

Proof. Since L is upper Hessenberg and A and B are strictly lower triangular, one 
immediately sees that LA — BL is lower triangular. Now L is also upper triangular, so 
that 

(Vi^LM" 1 )) L) = (f(LM~ 1 )l) 

and 

^(^(M- 1 ^)) \ =(lF(M- 1 L)) . 
This means that 

(la-bl)_ = (f(lm- 1 )l-lf(m- 1 l)} . 

From Lemma |3.1| we then obtain that (LA — BL)^ is the zero matrix. We have already 
mentioned that LA — BL is lower triangular, so we conclude that LA — BL is a diagonal 
matrix. 

Since M is lower triangular and A and B are strictly lower triangular, MA — BM 
is strictly lower triangular. Then we note that only the first of the subdiagonals of M 
has non-zero elements. This implies that below the first subdiagonal, 

MA - BM = (f (LM- 1 )) M — m(f (M~ x L) ) 
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agrees with 

F (LM- 1 ) M - MF (M~ X L) = O. 

The last equality follows from Lemma \5.l\ Thus the only non-zero elements of MA—BM 
are on the first sub diagonal. □ 



Remark 3.3 The generalized finite RLT (|3.2| ) can also be written as 

a n =(F(LM~ 1 )) -(f(M- 1 L)\ . 1< n < N, 

V / n.n— 1 V / 



n+l,n 

b„ = [ F (M _1 L) ) -(F(LM- 1 )) 



Kn< N -1. 



(3.5) 



Remark 3.4 If we take F(x) = 1/x in ( p.2|) we obtain the first representation of the 
finite RTL (|1.2|) . The other representation ( |1.3|) found in the literature, is obtained by 
taking F(x) = x. To see this, we note that from (|1.4j) we have 



M' 1 



( 1 

h 1 

b 2 bi bi 



0\ 





1 



V 



N-l 

n ^ 







bN-lbN-2 fcjv-i 1 



so that (LM = 6i(6i + i+a i+ i) and (M = 6 i (6 i _ 1 + a i ) for % = 1, 

So, if F(x) = x, the equations ( ^.5| ) are 

d n = b n -i{b n + an) - b n (b n -i + a n ) = a n (6 n _i - 6 n ), 



N-l. 



bn — b n (b n -x + a n ) — b n (b n+ i + a n+ i) — b n [a n — a n+ % + 6 n _ 



1 — "n+lj) 



with b = b N = 0. This corresponds to (|1.3|). 



5.^. T/ie evolution of the spectral data 

We now describe the time evolution of the spectral data introduced in Section |2|. The 
generalized RTL ( |3.2| ) can then be solved with the aid of this spectral transform which is 
a bijection between (M -1- )^" 1 and its image {0 < Ai < • • • < A^, Wj > 0, J2f=i w j = !}■ 
In particular we can start from positive initial matrix data a n (0) and b n (0) and follow 
the scheme in Figure |T] to find the matrix data at time t E R. 

Let L(t) and M(t) be the solutions of the generalized RTL (|3.2|) with 

A(t) = ~[f (M _1 (*)L(*))) and = -(F (L(t)M- x {t))) . (3.6) 

Let Ci and £2 be the unique solutions of the linear differential equations 

C 1 (t) = C x {t)B{t), A(0) = J, (3.7) 
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and 

C 2 (t) = -A(t)C 2 (t), £ 2 (0) = /. (3.8) 

Then £1 and £ 2 are lower triangular matrices with ones on the diagonal and therefore 
invertible. Furthermore 

j t {£iL£ 2 ) = &BLC2 + Cx{LA - BL)C 2 - C X LAC 2 = O 

and similarly |(£iM£ 2 ) = O. This means that 

C x {t)L{t)C 2 {t) = L(0) and C x (t)M(t)C 2 (t) = M(0). (3.9) 

Theorem 3.5 The eigenvalues \j of the generalized eigenvalue problem of the pair 
(L, M) are constants of motion, which means 

Xj = 0, 1 < j < N. (3.10) 
Proof. From fl3.9| ) it follows that for every z, 

C x (t)(zM(t) - L(t))C 2 (t) = zM(0) - L(0). 
Since det(£i(£)) = det(£ 2 (i)) = 1 we see that 

det(zM(t) - L(t)) = det(zM(0) - L(0)). 
This proves the theorem. □ 

To find the evolution of the spectral data Wj, we need the following lemma. For a 
matrix X we use X> to denote its upper triangular part. So X = X> + X-. 

Lemma 3.6 Let 1Z be the unique solution of the linear differential equations 

K(t) = -(F (L(t)Af _1 (*))) 7^(0, Tl(0) = I. (3.11) 

Then TZ is upper triangular and we have 

£ 1 (t)^(t)=e-* F ( L(0)A/ " 1(0) ), (3.12) 

which gives an LR- factorization of e _ * F ( L (°) M (°))_ 

Proof. Combining (|3.7f) , ( P-llQ , and (|3lj|) , we obtain 

j t (c x {t)n{t)) = £ x (t)B(t)K(t)-£ x (t)(F(L(t)M-\t))) > K(t) 
= - C x (t)F (L{t)M- l {t))K{t). 

From ( |3.9| ) we get that 

F {L{t)M-\t)) = F (£^ 1 (t)L(0)M- 1 (0)£ 1 (t)) 
= £r 1 (t)F(L(0)M~ 1 (0))£ 1 (t) 
(see 0, Definition 6.2.4]). Hence 

d (c 1 (t)K(t)) = -F (L(0)M-\0)) C x (t)K(t), 



dt\ 

which proves (|3~T2|) since C x (0)ft(0) = J. □ 
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The theorem below gives the evolution of the spectral data Wj. Moser obtained a 
similar expression for his spectral data in the case of the nonrelativistic Toda lattice 

Theorem 3.7 For the generalized RTL ( |3.2| ) with positive initial data a n (0) > 0, 
b n (fi) > 0, the spectral data Wj have the time evolution 

Wj (t) = -jp^ , l<j<N. (3.13) 

£ w k (0)e- tF ^ 

k=l 

Proof. Set D = diag(Ai, . . . , Xn), which is time independent. The matrix V(t) 
contains the right eigenvectors (with the first components equal to 1) of the generalized 
eigenvalue problem of the pair (L(t),M(t)). For t — this means that L(0)V(0) = 
M(0)V(0)D. Using Q) we then obtain 

L(t)£ 2 (t)V(0) = M(t)£ 2 (t)V(0)D. 

So the columns of the matrix £ 2 (t)V(0), which are £ 2 (t)vj(0), 1 < j < N, are right 
eigenvectors of the generalized eigenvalue problem of the pair (L(t), M(t)). Because 
£ 2 {t) is lower triangular with ones on the diagonal, and Vj(0) has first component equal 
to 1, we find that £ 2 Vj(0) has first component equal to 1. So we have 

Vj{t) =£ 2 (t)^-(0), l<j<N. (3.14) 

Similarly we get from U(0) T L(0) = DU(0) T M(0) and Q 

f/(0) T £i(t)L(t) = DU^f^^Mit). 

The jth column of the matrix £i(t) T U(0), which is £i(t) T Uj(0), is then a left eigenvector 
of the generalized eigenvalue problem of the pair (L(t), M(t)) with eigenvalue Xj. So 
it is equal to Uj(t) up to a constant factor. We now look at the first component of 
£i(t) T Uj(0). From Lemma |3.6| we get 

£ l (t)=e- tF ^ )M ' 1 ^)n-\t). 

Since lZ{t) is upper triangular, we obtain 

tf,(0) T A(f)gi = uM T e- tF{mM ~ H0)) (3.15) 

From the fact that Uj(0) is a left eigenvector of L(0)M _1 (0) with eigenvalue Xj, it follows 
that it is also a left eigenvector of e - tF ( L (°) M 1 (°)) with eigenvalue e~ tF( ~ Xj \ So 

^■(0) T e- tF ( L(0)M " 1(0) ) = e- iF ( A ^-(0) T 
Using this in (|3.15|) and noting that the first component of Uj(0) is 1, we see that 

^(0) T A(t)e!=e-^ {K- 1 ®)^. 
From this we conclude 

A(t) T %-(0) = e" iF ( A ^ {K-^t))^ uj(t), l<j<N. (3.16) 
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Combining Definition |2]3|, ([3.14Q , ( |3.16| ), and ( |3.9| ) we now obtain 
Wj(t)= fe(*) T M(t)%(t)) _1 

e tF(Ai) 



^(0) T A(i)M(i)£2(£K-(0) 



= {n-\t)) lfl v, j (o)e- tF M 

Since the sum to 1 we have 

fr-Wu-n • 

k=l 

which completes the proof of the theorem. □ 

Remark 3.8 Standard existence and uniqueness results from the theory of ordinary 
differential equations, give that the generalized RTL (|3.2|) with positive initial data 
a n (0) and b n (0) has a unique solution in a time interval containing t = 0, so that 
a n (t) > 0, b n (t) > 0. The expressions ( |3.10|) and (|3.13|) allow us to define spectral data 



for all t E R that remain in {0 < Ai < • ■ ■ < A^v, Wj > 0, J^jLi w j = ^ then follows 
that the system ( |3.2|) has a unique solution valid for all time, in which the matrix data 
are positive. 

3.3. Example 

As an example we take N = 5 and positive initial matrix data 



oi(0) 


= 3, 


fei(0) = 


1, 


a 2 (0) 


= 12, 


62(0) = 


6. 


a 3 (0) 


= 16, 


63(0) = 


11, 


o 4 (0) 


= 7, 


64(0) = 


5, 


a 5 (0) 


= 5. 







We obtain the corresponding spectral data, defined in Section pTT| , by solving the 
generalized eigenvalue problem of the pair (L(0),M(0)). The results are 



Ai 


= 1.9812757881, 


wi{0) = 


0.0097186754 


A 2 


= 2.6941860907, 


w 2 (0) = 


0.8409233539 


A 3 


= 6.6927423653, 


w 3 (0) = 


0.0757415291 


Ai 


= 13.8305993379, 


wt(Q) = 


0.0665694128 


A 5 


= 40.8011964181, 


w 5 (0) = 


0.0070470286. 



We consider the case that F(x) — 1/x. Knowing ( |3.10| ) and ( p. 13| ) we then compute the 



spectral data at several time values. For each of these values we construct the T-fraction 
( p.lOj ) as in the proof of Theorem [2.8| . This then gives the values of the matrix data. 



The results are shown in Figures || and ||. We observe that the b n tend to as t — > ±00. 
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Figure 2. Evolution of the matrix data a n (1 < n < 5). 



It is also clear from Figure ^]that the a n have limits for t — > ±oo. In the next section we 
will show that these limits are in fact Xq-u and X n as t — » +oo, or t — ► — oo, respectively. 



4. The sorting property of the generalized finite RTL 

In this section we suppose that F is strictly increasing and we investigate the behaviour 
of the matrix data of the generalized finite RTL when t — > ±oo. We will show that 



lim& n (t) = 0, l<n<JV-l, (4.1) 



and 



lim a n (t) = A n , lim a n (t) = Xn-u+i, 1 <n < N. (4.2) 

t— >+oo t^—oo 

We also establish the precise rates of convergence of ( f4.1| ) and ( [4.2[ ). For the SVD flow 
(and the related nonrelativistic Toda lattice) this was obtained in 0] . Note that in view 
of the sorting property in (|4.2|) the choice of spectral data Ai < ■ • ■ < Xn appears to be 
very natural. In the limit t — > ±oo, the diagonal entries a n (t) of the matrix L(t) tend 
to the eigenvalues and the eigenvalues are sorted in increasing order (as t — > +oo) or in 
decreasing order (as t — > — oo). 

The case when F is strictly decreasing can be studied in the same way as we will 
do for strictly increasing F. In this case, the b n (t) tend to 0, and the a n {t) tend to the 
eigenvalues in decreasing order as t — > +oo, and in increasing order as t — > — oo. 

Our proof of (|4.1|) and (|4.2|) is based on properties of the Laurent orthogonal 



polynomials P n . To indicate the time dependence we write P n (z,t). These polynomials 
are defined by the recurrence relation 

P n (z,t) = (z - a n (t))P n ^(z,t) -b n -i(t)zP n - 2 (z,t), (4.3) 
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Figure 3. Evolution of the matrix data b n (1 < n < A). 



with Pq = 1 and P_i = 0. For the polynomials P n (z,t) we find the following limits as 
t — > ±oo. 

Lemma 4.1 Assume F is strictly increasing. Then we have for 1 < n < N, 

n N 

lim P n (z,t) = T\(z-X k ) and lim P n (z,t) = TT (z - X k ). (4.4) 

t— »+oo t^—oo 

k=l k=N-n+l 



Proof. First we investigate what happens when t — > +oo. In Theorem |2.9| we 
showed that the polynomial P n (z,t) is orthogonal with respect to the discrete measure 
Y^j=i~0^^j- ft * s well-known that monic orthogonal polynomials minimize the 
weighted L 2 norm, see [O, Section 3.1]. In our case this means that the minimum 
of 

£(^))^ («) 

j = l 3 

taken over all monic polynomials q n of degree n, is attained by P n (z,t). Take 
q n (z) = Ylk=i( z ~ then (since Xj > and Wj(t) > 0) we obtain for 1 < i < N, 



3=1 3 



< > iaJX^^ 



JV 



j=n+l k=l 3 



The relativistic Toda lattice and Laurent orthogonal polynomials 20 
From Theorem |3.7| we get for the ratio of the spectral data Wj(t) that 



^ A-'{ FM - F ^, l<j,k<N. (4.7) 
Wk{t) w k {0) 

Since F is strictly increasing, we obtain from 

lim ^4-7 = 0, 1 < k < j < N. (4.8) 

t^+oo Wk{t) 

So (|4.6|) gives for every 1 < i < N, 

(p n (A,,t)) 2 ^< (n(A n+1 -A fe )) 2 ^±i^(l+o(l)), ast^+oo. (4.9) 

i j, = i n+1 

For 1 < % < n, we get from ([4. 81) that Wn+ }}^ tends to as t — > +oo. So from ((4.91) we 



„, ..^ g,^ -vw w ^ 

conclude that 

lim P n {Xi,t) = 0, 1 < i < n. 

t— >+oo 

Then the first limit in ( (4.4| ) follows, since P n (z,t) is a monic polynomial of degree n for 
every t . Similarly we can prove the second limit when t — > — oo, by taking 

N 

i=AT-n+l 

in (JO). □ 



Starting from the recurrence relation (|4.3j ), the limits of the polynomials P n (z,t) 
as t — > ±oo give rise to the limits of the matrix data ( |4.1D and (|4.2j ). However, we also 
want to establish the precise rate of convergence. We then need the following lemma. 

Lemma 4.2 Let < n < N - 1. Then for 1 < I < n + 1, 

n / n+1 \ /x 

J] &*(*) = P n {\ u t) A t ) (l + o(l)j , as t - +oo, (4.10) 

fc=l \i = l,2^Z / * 

and /or N — n < I < N , 

f[h(t) = P n (X h t) ( f[ (A,-Ai))^ (l + o(l)), ast^-oo. (4.11) 

fc=l \i=N-n,i^l / * 



Proof. From Theorem |2.9| it follows that for every monic polynomial Q n of degree n, 
we have 

N , s n 

^P„(A„t)Q„(A,)^ = n^W> 
j=i i fc=i 

where < n < N — 1. If we take 1 < I < n + 1 and 

n+1 

Qn(»= n ^" A *)' 

i=l,i^ 
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f[b k (t) = P n (X h t) ( f[ (X l - Ai) ) 

k=l 
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A? 



a; 



+ j2 n ^-^ 



j=rc+2 



From Lemma 14.11 we see that 



lim P n (\j,t) ^ 0, n + 1 < j < N. 



t— >+oo 



(4.12) 



(4.13) 



Taking I = n + 1 in ( 14.121) and letting t -> +oo, we get with ( |4.13|) and O that 

f[b k (t) = P4X n+u t)(f[(X l -X i ))^^ (l + o(l)), ast^+oc, (4.14) 

k=l \i=l J l 

n 

which is ( f4.1(J| ) with I = n + 1. So, for large t, Yl bk(t) behaves like w n+ \{t) up to a 

fe=i 

constant. Knowing this, we get (|4. 10|) for 1 < / < n from (|4.12|) , (|4.13| ) and (|4~8|) . The 
proof of (|4.11| ) is similar. □ 



From Lemma we can easily establish the rates of convergence of P n (Xi,t), 
1 < I < n, as t — > +oo, and of P n (A/, t), N — n + 1 < I < N, as t — > — oo. 

Now we look at the matrix data a n (t) and b n (t). Lemmas £T and L2 will help us 
to find their behaviour when t — > ±oo. 



Theorem 4.3 For the generalized finite RTL (|3.2| ) with F strictly increasing we have 
fori <n < N- 1, 

2 



b n (t) 



w n+l (0) (A n ) 



n-l 



w n (0) (A n+ i) 
as t —* +oo, and 



A 



EI (An+l — Aj 

i=l 

n-1 

^ I! i^n - Aj) j 

( 



-t[ F(X n+1 )-F(X n ) 



l + o(L 



(4.15) 



b n {t) 



W N _ n (0) (Aat„ 



in-1 



-n+1, 



Wat_„ + i(0) (Ajv- 



iV \_ 

n {Xn-u — Ai) 

i=N~n+l 
N 

, II (AjV-n+l - A^ 
\i=AT-n+2 



x e 



t \ F(X N - n+1 )-F(X N - n ) 



l + o(l) , (4.16) 



as t ^ — oo. 



Proof. We only give the proof for the case when t +oo. The proof of the case 
t — > — oo is similar. Let < n < N — 1. From Lemma PO we have that 



lim P n (X n+ i,t) = TT(A n+ i - Ai 



t— >+oo 
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So taking I — n + 1 in ( |4.10| ) we obtain 

\ 2 

\ in , J+\ / \ 

as t -> +oo. (4.17) 



n 6 *(*)=fn^-^) ^7 (i+o(i) 



Since 



W n+1 (t) _ W n+1 (0) -t(F(\ n+1 )-F(\ n )j 

w n (t) w n (0) 

( |4.15| ) follows immediately from (|4.17|) . □ 

This theorem implies that the matrix data b n (t) tend to as t — > ±oo at an 
exponential rate. We now use the recurrence relation for the polynomials P n (z,t) to 
compute the behaviour of a n (t) as t —>■ ±oo. 

Theorem 4.4 Let 1 < n < N. For the generalized finite RTL ( |3.2| ) with F strictly 
increasing we have 

A n - a n {t) = \ b n (t) (l + o(l)) - - An - b n ^(t) (l + o(l)) , (4.18) 
as t — > +oo, and 



Ajv-n+i - a n (t) = - — ^ N n+1 b n {t) (l + o(l)) 



A 



iV-n+1 



Av-n+2 — ^N-n+l 

as t — > — oo. Recall that b^ = 0, 6q = 0. 



6 n _ 1 (t)(l+o(l)), (4.19) 



Proof. We will only prove ( |4.18|) and in the proof we assume 2 < n < N — 1. The 
recurrence relation fl4.3|) for the polynomials P n (z,t) gives 

An - a n (0 = A ^ (A " ;t) b n (t) + K P ;-f^ b n -i(t). (4.20) 
From Lemma |4.1j we obtain 

l^ -T-Vfl+oW), ast^+oc. (4.21) 



From (|4.10|) we have 



n+l x 1 



Pn l Aw,t) tonft) = fl (An-A,)) (l + o(l)Y ast->+oo. 
K II h(t) \<=i,¥n / 



fc=l 

and 



-fft-i(A n , t) 
w ' 

(a.)™- 1 n 



(t) = ( H(K- A,) j (l + o(l)), ast^+oo. 



Taking the ratio of the last two expressions, we obtain 
-fn(A n ,t) 1 



A n b n (t)P n -i(\ n ,t) X n — A n+ i 



(l + o(l)), ast^+oo. (4.22) 
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Combining (|O0|) , (|OI|) and ( ggg) gives (fO|). 

The cases n = 1 and n = N are simpler, because then either the last or the first 
term on the right-hand side of ( 4.20| ) equals 0. The proof for the case t — > — oo is similar. 



□ 



Since we know from Theorem [13] that the b n (t) tend to as t — > ±oo, Theorem [O 
proves the limits ( |4.2j ). Because the b n (t) are exponentially small as t — > ±oo, we also 
see that a n (t) approaches its limit at an exponential rate. 



Remark 4.5 We return to the Newtonian form (.LI) of the finite RTL. In particular 
we look at the behaviour of the differences of the q n as t — > ±oo. In |], |2| it was shown 
that the change of variables 



h(q, 



n-l 



,Pr, 



h(q n - <?n+l) 

2 exp(q n - q n+1 + p„)/i(g n _i 



Qn) 



l<n<N, 
1< n< N- 1, 



(4.23) 



where h(x) = vl + e 2 e x , leads to ( |1.3| ). This system is of the form ( p.2| ) with F(a;) = a; 
(see Remark |3.4| ). Here the initial conditions a„(0) and 6 n (0) can be computed from the 
q n (0) and p n (0) (or g n (0) and g n (0)). From ( |4.23| ) we have b n = e 2 e~ ( - qn+1 ~ Qn ' , a n , so that 

... , ,'e 2 a„,(±oo) 

QVm-iW - ^nW = log 



b n (t) 



o(l), 



as t — > ±oo. 



The asymptotic results ( |4.1| ) and ( |4.2j ) for the matrix data then show that the differences 
between the q n , namely q n +i — q n , tend to +oo for t — >• ±oo. From Theorem [4.3| we also 
know that the b n decay at an exponential rate as t —>■ ±oo. So the g n+ i — q n have linear 
behaviour as t — >• ±oo, where the intercept and the slope can be established explicitly. 
In particular we get 

?«(*) =*(A n+ i-A n ) + 0(l), 



?n+l(*) 

and 

9rH-l(*) 

where the values Ai < 



as i 



-oo, 



9n(*) = (Aiv-„ + i - Aat_„) + 0(1), as t -> -oo, 

• • < Aat are obtained from the initial conditions a n (0) and 6 n (0). 
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